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Abstract 

We consider a five dimensional vacuum cosmology with Bianchi 
type-IX spatial geometry and an extra non-compact coordinate. Find- 
ing a new class of solutions, we examine and rule out the possibility 
of deterministic chaos. We interpret this result within the context of 
induced matter theory. 

1 Introduction 

The Mixmaster universe, an anisotropic, homogeneous cosmology with Bianchi 
type-IX spatial geometry, has long been a subject of interest. Belinskii, Kha- 
latnikov and Lifshitz |l| analyzed its dynamics in their studies of the general 
approach to the initial singularity. Misner proposed a period of Mixmaster 
dynamics as the possible solution to the horizon problem. 

One interesting feature of these models is deterministic chaotic behav- 
ior. This manifests itself in a series of "Kasner epochs" in which two of 
the models' scale factors oscillate, approximated by Kasner solutions linked 
asymptotically, while the third scale factor monotonically expands. After a 
period of time, a so-called "era," the expanding scale factor switches places 
with one of the oscillating factors - the second, for instance. The dynamics 
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then continues with the first and third factors now oscillating, and the second 
monotonically expanding. This progression from one era to the next occurs 
indefinitely, with the number of epochs in each successive era appearing as 
random as a dice toss. By statistical measures, this represents deterministic 
chaos, as Barrow || has shown and many others have confirmed ||, [5|, |, [7]]. 

Chaos is found in the standard 4D vacuum Bianchi type-IX universe, as 
well as in the vacuum Bianchi type- VIII. However, for all but a few cases, 
once matter and energy fields are included (and these are too large to ignore; 
i.e., not small perturbations), chaotic behavior is suppressed [§, |§. Moreover, 
as Barrow and Stein-Schabes [M, as well as Furusawa |TT| and others |L2 



|T3| , |H|1 (working independently using various methods) have demonstrated, 
chaos is absent even in vacuum cases for models of five dimensions or higher. 
The absence of higher dimensional chaos includes brane world models, as 
Coley[T3|] has recently determined. 

Might there be a connection between the suppression of chaos in matter- 
filled 4D models, and the absence of chaos in vacuum 5D (or higher) models? 
This linkage can be examined by means of induced matter theory [|16|] , which 
associates the additional geometrical terms arising from the extra scale factor 
and derivatives of the fifth coordinate in the five dimensional vacuum Einstein 
equations with the matter-energy components of four dimensional theory. In 
other words, 4D matter stems from 5D geometry. 

To investigate this question, we examine the properties of Bianchi type-IX 
models with the addition of an extended fifth coordinate. By calculating the 
density and pressure of the induced matter, we demonstrate that this model 
cannot be regarded as empty. Therefore, its properties are qualitatively 
different from the standard vacuum, non-Kaluza-Klein 4D result, and do not 
manifest chaos for this reason. 

In fact, according to induced matter theory, there is no pure 4D vacuum, 
as can be seen from the canonical form due to Mashhoon, Liu and Wesson 
[IT7]. In the 5D Bianchi type-IX case, this precludes chaos, as we demonstrate. 



2 5D Mixmaster Dynamics 

We consider the five dimensional metric: 

ds 2 = e u dt 2 - gij Ju j - e^dl 2 (1) 
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as: 



where the 3D spatial part of the metric can be expressed in diagonal form 

g ij = diag(e a ,el 3 ,ei) (2) 



The time coordinate t and the three spatial coordinates x, y and z have 
been supplemented with a fifth coordinate I. We assume that the metric 
coefficients //, u, a, f3 and 7 each depend, in general, on both t and /. 

The one-forms iv % have the relationship: 

du l = ho i jk oo ] uj k (3) 

where the non-zero structure constants are C 1 2 3 = — C ,1 32 = 1, and cyclic 
permutations. 

The 5D Einstein tensor for the Bianchi type-IX case can be written in 
the form: 

111111 

e "( 2 2 2 « 

-/x a - -/x /3 - -ft 7 + -a /3 +-«7 + -/3 7 ) + 

e" e 7 1 1 1 

4 e /3 e 7 + Ae a e"< + Ae a e< 3 ~ 2e^~ 2e^~ 2e< ^ 

G° 4 = e~ v { 2a-* + 2/3"* + 2 7 '* + a'a* + + 7-7* - 

a u* - (3-u* - r»* - «>■ - /?> - 7» ( 5 ) 
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—v I ^ .. ^ .. ^ .. ^ .2 1 .2 1 .2 
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2 v 2 2 2 4 4 4 
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2 2 2 4 4 4 

-a p + -a z/ + -/? z/ - -a /i - -By, — -a v ) — 
4 4 4 4 r 4 r 4 r/ 

e Q 3e 7 1 1 1 

4e /3 eT 4e Q e^ 4e a e' 3 2e a 2e< i 2e< 

Gl = e -( -ia-^-ir-i^-i^-ir'- 

-arp' - T a '7' ~ -7 Pi' + t"'^' + -,P' V ' + tT'^'J + 
4 4 4 4 4 4 
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(9) 



4e^e^ 4e a eT 4e a e^ 2e a 2c/ 3 2eT 

where we use overdots to represent partial derivatives with respect to t, 
and asterisks to represent partial derivatives with respect to I. 

Following the standard procedure, we collect each of the terms in the 5D 
vacuum Einstein equations, C„ = 0, dependent on either ji or on derivatives 
with respect to /, place them on the right-hand side and identify them as the 
induced matter density and pressure components: 
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-// a - -n (3 - -n 7 + -a /3 + -a 7 + 7 XlO) 
1 1 , 1„ 1 1 

2 P 2 7 2 4 P 4 7 4 

1 1 2 1 1 In 

8ttP 2 = e ( -/i- + -fi- + -a>- + ~rfi m - ~^ i V ') + 

e-"( --a** - -7** - V* - -a* 2 - - V 2 - 
v 2 2 ' 2 4 4 ' 4 

-a 7 - -a 1/ - -7 ^ + /i + ^7 A* + ^ v )( 12 ) 

1 1 2 1 L 1 

8ttP 3 = e ( -fi" + + -a-p + -pp - -^i v ) + 

111111 

e -N — a ** _ -R** _ _ _ a * 2 _ _/?* 2 _ _^* 2 _ 

v 2 2 P 2 4 4 P 4 

- |oV - ±/JV + I«V + + ^V*)(13) 

The Einstein equations now can be written in 4D form as: 

e-"( --«•/?■ - -a'7 - + 

e Q e^ 3 1 1 J_ _ 

4 e /3 e 7 4e a eT 4e a e^ 2e a 2e^ 2e^ ^ ~~ 

«-( - \r - Y - - \?r + \t* + xrr) + 

3e Q e 1 1 1 1 

+ ^"^-^-8^=0 (15) 



4e^eT 4e Q eT 4e a e /3 2e a 2e^ 2eT 
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1 1 1 2 1 2 1 1 1 

--a" 7" a' 7' a'7' H — ocv H — 7'f 

2 2' 4 4' 4 ' 4 4 ' ; 

3e^ e 7 1 1 1 



4e/3 e 7 + 4e a e^ 4e a e£ 2e a + 2e^ 2e^ ^ ° ^ 

e~"( —a" - -8" - -a- 2 - -0- 2 - \oc& + \ocv + -B'v) - 
v 2 2 4 4 4 4 4 ; 

e a e 13 3e 7 1 1 1 .. 

4e%^ ~ Ae a e"i + 4e a eP ~ 2e« ~ 2e? + 2e^ ~ ^ 3 ~ ' ^ 

e _ "( 2a-* + 2/?-* + 2 7 -* + era* + 8'8* + ri* ~ 

a v* - Q-v* - rv* - a* ^ - (3*ll- - 7>') = (18) 

-vl 1 .. 1 a.. ^ .. 1 .2 1 o.2 1 .2 

e — a p 7 a p 7 — 

v 2 2 M 2 ' 4 4 P 4 ' 

1 a 1 1, 1 L 1 

-a-p - -a'r - -fJ'T + 2 a ' u ' + \P' V ' + 4^' z/ ') + 

111111 

e **( -a /5 + -a 7 + -/5 7 + -a i/ + -/3 i/ + -7 z/ ) + 

e " 67 1 1 1 - fl9) 

4 e /%7 + 4 e «eT + 4e a e^ ~ 2e" ~ 2e^ ~ 2e^ ~ ' ' 

We solve equations ( |i~4] - |i~9| ) through the method of separation of variables. 
We write: 



e a = [a(t)a(l)] 2 (20) 

e P = W)b{l)f (21) 

e 7 = [c(t)c(l)} 2 (22) 

= \d(t)d(l)} 2 (23) 

e" = [f(t)f(l)] 2 (24) 

where we have explicitly indicated dependence on either t or I exclusively. 
We then make the following simplifying assumptions: 

/ = abed (25) 
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a = 


e sil 


b = 


e S2l 


c = 


e szl 


d = 


e S4l 


f = 


e s l 



(26) 
(27) 
(28) 
(29) 
(30) 

where the Sj are all constants. 

To solve, we impose the conditions: 

si = (31) 

s = s 2 + s 3 (32) 

Equations (|15|-|T7j) then reduce to: 

(^)- = I((6V S2i -cV^) 2 -« 4 )^ 2 (33) 
(t) = I(( c V S3 '-a 2 ) 2 -6V S2 )rf 2 + 

(2s 2 2 + 2s 2 s 3 - s 2 s 4 )a 2 b 2 c 2 e 2 ^ 2+s3 -^ 1 (34) 
(-)' = l -{{a 2 -b 2 e 2s *) 2 -c i e^)d 2 + 

{2s 3 2 + 2s 2 s 3 - s 3 s 4 )a 2 b 2 c 2 e 2( - s2+s3 - s4 » (35) 

Following the procedure of Belinskii, Khalatnikov and Lifshitz, we exam- 
ine this system's evolution backward in time. If one could neglect all the 
terms on the RHS of equations ( |33| - |35l) then one would obtain a 5D general- 
ized Kasner solution, similar to that found by one of us in an earlier paper 



lfil, with: 



a ~ e Apit (36) 



b ~ e Ap2t (37) 
c ~ e Ap3t (38) 

where A and Pi (the generalized Kasner exponents) are constants. 
The relationship amongst the generalized Kasner exponents is: 



i 



EP* = 1 ( 39 ) 

i=i 

E^ 2 = 1 (40) 
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Suppose, as these constraints allow, one of the three spatial exponents is 
negative-pi, for instance-and the other two, p 2 and p 3 , are positive. In that 
case a would increase as the system evolves backward in time, and b and c 



would decrease. Soon the RHS of equations (|33| - p5D would become dominated 
by the terms containing a 4 , and they would reduce to the approximate form: 



a 

,b_ 

"b 

,c 



1 4 j2 

--a a 
2 



4 J2 



-ad 



4 j2 



-ad 



y c' 2 

In addition, equation (|19D can be written as: 

(±) = 2(s 2 2 + s 3 2 + s 2 s 3 - s 2 s 4 - s 3 s 4 )e 2 ^-^ l a 2 b 2 c 2 



which yields: 

d ■■ 

where A is a constant. 
Substituting (g5|) into (gTJ-| 
yields: 



S2 2 + S 2 S 3 + g 3 2 
S2 + S 3 



(41) 
(42) 
(43) 

(44) 

(45) 
(46) 



and then solving along with equation (|I~4|) 



A(2pi+p 4 ) _ e -l/2A P4 t 



b 

c 



\ cosh (A (2 Pl +p 4 )t) 

b ^COsh (A (2 Pl + p 4 ) t ) e -l/2A*(-2 P 2+2 Pl+P4 ) 
C ^COsh (A (2 p X + p 4 ) t ) e -l/2A t (-2 P3 +2 Pl + P 4) 



(47) 

(48) 
(49) 



where bo and Co are constants. The choice of remaining constants is such 
that for large t the solution asymptotically matches the generalized Kasner 
behavior described by (PE1S8|). 
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Taken into account conditions (|3T| , ^ and [IS), as well as equation (|T8|), 



we can now express the full solution as: 





= A(2pi + p 4 )&o 2 c 2 cosh (A(2pi + p 4 )t) 




e (A(-4pi+2 P 2+2p3-P4)t+28 3 (l-'e)0 


e a 


A (2p!+p 4 )e- Ap4 * 




cosh (A (2pi + p 4 )t) 




= 6 2 cosh (A (2 Pl + p 4 ) t) e ( 2 A*w-2At Pl 


e 7 


= c 2 cosh (A (2 Pl + p 4 ) t) e ( 2A ^- 2A ^ 




__ e (2Ap 4 t+4s 3 


where k = 


|±|V5. 



(50) 

(51) 

(52) 
(53) 
(54) 



Note that for S3 = 0, dependency on Z vanishes, and this solution reduces 
to that of Barrow and Stein-Schabes ||10| . Otherwise, for non-vanishing S3, 
we obtain the additional constraint: 

Pi -piK -p 2 +p 3 n = (55) 

Let's now examine the epoch to epoch behavior of this model. For t — > 00 
its spatial scale factors asymptotically match the generalized Kasner solution: 



a = 2Apit (56) 
(3 = 2Ap 2 t-2Ks 3 l (57) 
7 = 2Ap 3 t + 2s 3 / (58) 

(59) 



Then, as t runs backward, its scale factors evolve according to ([51] - [53]) . 
Taking t — > —00 (£ in this coordinate system ranges from -00 to 00), its spatial 
scale factors asymptotically match another generalized Kasner solution: 



a' = 2A' Pl 't (60) 
f3' = 2A'p 2 't - 2ks 3 1 (61) 
V = 2A'p3^ + 2s 3 / (62) 

(63) 
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where the old and new epochs are connected by the relationships: 



AW 
AW 
A'W 
AW 



-A(p x + p 4 ) 
A(p 2 - 2pi - p A ) 
A(ps - 2pi - pa) 
Apa 



(64) 
(65) 
(66) 
(67) 



These epoch-epoch transitions continue until pi,p 2 and p 3 are all positive. 
Such a circumstance is permitted by conditions fl39[ , |40| and ^) for a finite 
range of generalised Kasner exponents. When that happens, the RHS of 
equations ( |33| - |35l) vanish, and no more transitions take place. The universe 
continues to evolve indefinitely with the same set of exponents. Because the 
number of transitions is finite, this solution is manifestly not chaotic. 

3 Induced Matter Properties 

Let's now examine the density and pressure of the induced matter associated 
with this solution. Substituting into eqns. (|Il3|-|i"3|) , we obtain: 





p 4 2 v/-20p 4 ^ + 10p4 + 10 



(68) 




p 4 2 V-2<W + 10p 4 + 10 



(69) 




p 4 2 V-20j9 4 2 + 10p 4 + 10 



(70) 




p 4 2 V-20j9 4 2 + 10p 4 + 10 



(71) 



where: 
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9i = fihh (cosh (l/3 A (p 4 + 2 + ^/-4p 4 2 + 2^4 + 2^) *)) (72) 
<? 2 = fj 2 f 3 (sinh (l/3 A (p 4 + 2 + + 2p 4 + 2y/E^ tjj (73) 



with: 



3 A tp 4 +6 A t (-1/3 P4+1/3+1/6 v /-20p 4 2 +10p4+10^ -2 A t+s 3 Z-s 3 J a/5 



/i = e """V — v " ' -V ' (74) 

, = ^ 

2& V (l/3 p 4 + 2/3 + 1/3 V-20V + 10p 4 + 10) 

/ 3 = sech 2 (a (l/3p 4 + 2/3 + 1/3 ^/-20p 4 2 + 10p 4 + 1(d) (76) 

Though the pressure is clearly oscillatory and anisotropic, we can define 
an effective function: 

^// = ^E^ (77) 

This yields the "equation of state" : 

p = 3P eff (78) 

which is that of a hot photon gas, expanding anisotropically, as may have 
been the case in the very early universe. 

4 Conclusions 

We have examined the behavior of a 5D universe with spatial slices that 
possess Bianchi type-IX geometries and a non-compact fifth coordinate. We 
have shown that this solution oscillates between generalized Kasner epochs, 
but not indefinitely, and therefore does not possess deterministic chaos. An- 
alyzing this solution within the context of induced matter theory, we have 
found an effective equation of state, resembling an anisotropically expanding 
photon gas. In this manner, induced matter theory has provided a possible 
explanation for the lack of chaos in vacuum 5D Bianchi type-IX models: as 
expressed in 4D form, the induced matter acts to suppress transitions. 
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